

clear
clear matrix
clear mata
set maxvar 10000

set more off

cd "/Users/uctqa20/Dropbox/Land&UrbanizationProject/Analysis/Replication Material"

use "harris_honig_jop_replication_data.dta"

*replace as missing any estimates that are based on less than 10 respondents
replace prop_titlers = . if sa_totpop < 10

collapse (mean) prop_titlers mal_border, by(sqkm)

lab var prop_titlers "Proportion in Community with Title"

*** Figure B1 (a) ***
hist prop_titlers if mal_border == 1, bcolor(gray) graphregion(color(white)) ///
xscale(range(0 .5)) xlabel(0(.1).5) yscale(range(0 25)) ylabel(0(5)25)

*** Figure B1 (b) ***
hist prop_titlers if mal_border == 0, bcolor(gray) graphregion(color(white)) ///
xscale(range(0 .5)) xlabel(0(.1).5) yscale(range(0 25)) ylabel(0(5)25)
